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ABSTRACT 

Bivariate luminosity functions (LFs) are computed for galaxies in the New York Value-Added 
Galaxy Catalogue, based on the Sloan Digital Sky Survey Data Release 4. The galaxy prop- 
erties investigated are the morphological type, inverse concentration index, Sersic index, ab- 
solute effective surface brightness, reference frame colours, absolute radius, eClass spectral 
type, stellar mass and galaxy environment. The morphological sample is flux-limited to galax- 
ies with r < 15.9 and consists of 37,047 classifications to an RMS accuracy of ± half a class 
in the sequence E, SO, Sa, Sb, Sc, Sd, Im. These were assigned by an artificial neural network, 
based on a training set of 645 eyeball classifications. The other samples use r < 17.77 with a 
median redshift of z 0.08, and a limiting redshift of z < 0.15 to minimize the effects of evo- 
lution. Other cuts, for example in axis ratio, are made to minimize biases. A wealth of detail is 
seen, with clear variations between the LFs according to absolute magnitude and the second 
parameter. They are consistent with an early type, bright, concentrated, red population and a 
late type, faint, less concentrated, blue, star forming population. This bimodality suggests two 
major underlying physical processes, which in agreement with previous authors we hypothe- 
size to be merger and accretion, associated with the properties of bulges and discs respectively. 
The bivariate luminosity-surface brightness distribution is fit with the Choloniewski function 
(a Schechter function in absolute magnitude and Gaussian in surface brightness). The fit is 
found to be poor, as might be expected if there are two underlying processes. 

Key words: cosmology: observations - methods: data analysis - methods: statistical - galax- 
ies: fundamental parameters - galaxies: luminosity function - galaxies: statistics 



1 INTRODUCTION 

The galaxy luminosity function (LF) is well-known as a funda- 
mental measurement of the properties of galaxies. It constrains 
theories of the ir formation and ev olution, such as biased galaxy 
formation (e.g. Benson et al. 2000) and is needed for many other 
mea surements. Examp les are the luminosity density in the universe 
(e.B. ICross et alJlOOJ). the selection func tion in magnitude-limited 
galaxy surveys (e.g. iNorberg et ai]|2002h . the distances to various 
types of object, the number of absorbing objects subdivided by red- 
shift and deprojecting the angular correlation function from two di- 
mensions to three (e.g. iLimber 1953.:,Binggeli etal, 1988 ). 

The general or universal luminosity function, i.e. that for all 
galaxies, is defined as the comoving number density of galaxies 
from luminosity L to L + dL 



dN = cj>{L) dL dV. 



(1) 



Measu rements of the universal form began with iHubbid 
who claimed a Gaussian form. Various studies fol- 
lowed, which included data on dwarf galaxies in clusters (e.g. Abd| 
1962) and the study of iHolmberd (l974|) of galaxies in the field. 
Both of these argued for a power-law faint-end slope to the overall 
LF, as opposed to a Gaussian form. 

In the mid 1970s it was suggested ^ Schechte3ll97d) that the 
optical LF can be approximated by the function 

<P{L) dL ^ c^* ( ^) e^p(-^)d(^) , (2) 
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where 0* is the normalisation, L* is the characteristic luminosity 
above which the function has an exponential cutoff and a is the 
value of the power-law slope below L* . The Schechter function, as 
it is now known, has been used for most subsequent characterisa- 
tions of the universal LF. The function i s motivated by self-simila r 
gravitational condensation of structures jPress & SchechteJl974) . 
except that a is not a fixed parameter. 

Until recently the study of the LF, particularly for field galax- 



© 2006 RAS 



2 A^. M. Ball et al. 



ies as opposed to those in clusters, has been hampered by the lack 
of large samples. Few galaxies had redshifts and the photometry 
was from photographic plates, resulting in a large variation in the 
measured LF between different surveys. In the past decade data 
from large red shift surveys, such as the 2dF Galaxy Redshift Sur- 
vey (2dFGRS. IColless et all200lLl2003^ and the Sloan Digital Sky 
Survey (SPSS. lYork et alj|2000h have become available, greafly 
improving the situation, although the LFs are still not in perfect 
agreement — for example. Xiske et al. I .2003.) gives revisions to the 
normalisations of the LF in various surveys based on the deeper 
wide-field imaging obtained by the Millennium Galaxy Catalogue 
(MGC) , which they also describe, and recent comparisons of LFs 
such as lDrivej E0041 conclude that the variations between surveys 
are dominated by the systematic errors, particularly at the faint end. 
The differences are most likely due to surface brightness selection 
effects. 

Whilst the universal LF is reasonably well constrained, it is 
well-known that the LF varies as a function of intrinsic galaxy prop- 
erties. Thus the universal LF needs to be augmented by a descrip- 
tion of this variation. As above, the lack of data until recently has 
hampered such studies, but with the large sample sizes (2dFGRS 
and SDSS) and five band CCD photometry (SDSS) now available, 
detailed divisions of the galaxies into samples which are still statis- 
tically significant is now possible. One can also apply fairly strict 
sample cuts to minimize biases and still retain large samples. 

Previous studies of non-universal LFs have included 1) those 
focusing on specific environments, particularly clusters as the 
galaxies are all at approximately the same distance and in a nar- 
row field of view, 2) those focussing on specific types of galaxies, 
the large number of studies which divide the general LF according 
to some criterion and fit functions to each bin, and 3) those which 
study the LF bivariate with another parameter, either binned or as 
an anal ytical function. Recent examples, so me of which are bivari- 
ate, are |PoDesso et al.' EoO?) for clusters, 'Mercuric et al.' ('2006') 
for su perclusters, Hovle et al. (2005) for voids, de Jong & Lacey 
(2000) and de LaDoarent et al. ( 200 4i) for spirals, a nd iReda et alJ 
t2004.) . .Stocke et aL(.2004..) and .de Jong et alJj20o4 for ellipticals. 

Various theories and simulations exist pertaining to the origin 
of the LF, tied in with the physics of galaxy formation. Examples 
are^enson et al. (2003) who begin with the mass function of dark 
matter haloes in the ACDM cosmology and add gas cooling, pho- 
toionization, feedback (e.g. from supemovae), galaxy merging and 
thermal conduction. Mo et al. ( 2004) assume that the segregation of 
the galaxy population by environment is due to that of halo prop- 
erties and use the h alo occupation distribution to g ive predicted 
LFs by environment. ICoorav & Milosavlievid i2005h reconstructs 
the Schechter form of the LF using empirical relations between the 
central galaxy luminosity and halo mass, and the total galaxy lumi- 
nosity and halo mass. 

Further details of the LF are in v arious reviews, for example 
iBinggeli. Sandage & Tammannlj 19881) and subdivided by morpho- 
logical type in lde^apDMenn2003l) ~^ 

In this paper the bivariate LF is studied in bins us- 
ing the nonparametric stepwise m aximum likelihood method of 
lEfstathiou. Ellis & Peterson il988l) and two-dimensional analyt- 
ical functions are fitte d to the resu l ts. Pr evious studies that 
have done this include 'Choloniewski' f 19851) subdivided by ra- 
dius, So dre & Lahav (1993) subdivided by galaxy diameter and 
ICross et al . (2001), Cross & Driver (2002) and Driver et al. (2005) 
subdivided by absolute effective surface brightness. The latter 
qua ntity is also known as the bivariate brightness distribution. 

ball et ai] i2004h showed that morphological types can be re- 



liably assigned to galaxies in the SDSS using artificial neural net- 
works (ANNs) provided a representative training set is available. 
Here types are assigned to the resolution E, SO, Sa, Sb, Sc, Sd and 
Im for 37,047 galaxies. In Ball et al. (2006, in preparation, hereafter 
B06) we use these types to study the variation of galaxy morphol- 
ogy with environment and compare the results to those for colour. 

The bivariate LFs are computed for ga laxies in the New Yor k 
Value-Added Galaxy Catalogue (VAGC, Bla nton et all t2005a), 
based on the Sloan Digital Sky Survey Data Release 4 (DR4, 
lAdelman-McCarthv et al. 2006, http://www.sdss.org/dr4i for a 
range of galaxy properties, including the morphological type de- 
scribed, inverse concentration index, Sersic index, absolute ef- 
fective surface brightness, eClass spectral type, reference frame 
colours, absolute 90% radius, stellar mass and galaxy environment. 
The galaxy samples are flux limited to r < 17.77, with the excep- 
tion of the morphological type, which is limited to r < 15.9. 

Here the LF is only considered in the wavebands of the 
SDSS, but there are m any studies in other wavebands, includ- 



ing t he ultraviolet (e.g. iBaldrv et ai]l2005t iBudayari et"an 20051: 
Wvd er et alj l2005h. near-infrared K band (e.g. [ Lovedavl boool 
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Cole et al.'200l':'Kochanek et al.'200lL lEke et alj2004 'jones et al 

2006), radio (e.g. Sadler etal. 2002) and x-ray (e.g. Mivaii et 
2000; Ranalfietal. 2005). 

.Blanton et al.. (.2003&) andlLovedavl j2004 show that the LFs 
in the SDSS DRI iAbazaiian et aj2003h to 2 < 0.3 are consistent 
with evolution in the galaxy population over this redshift range. 
Lovedav (2004) shows that even to a redshift of 0.15 evolution in 
the number density occurs, but that it is only marginally significant 
in DRI. Here, because we are measuring and fitting many different 
bivariate functions, we do not attempt to account for evolution, but 
restrict the sample to z ^ 0.15 to minimize its effect. 

Throughout, the standard spatial geometry is assumed, with 
Euclidean space, S^mattcr = 0.3 and Q.a ~ 0.7. For compat- 
ibility with previous studies, the dimensionless Hubble constant, 
h = JJo/lOO krns"^ Mpc"\ is set to 1. 



2 DATA 

The SDSS is a project to map n steradians of the northern galac- 
tic cap in five bands (u, g, r, i and z) from 3,500-8,900 A. This 
will provide photometry for of order 5 x 10^ galaxies. A multi- 
fibre spectrograph will provide redshifts and spectra for approxi- 
mately 10^ of these. A technical summary of the survey is given i n 
iYorketalit200a) . The telescope is d escribed iaGunn et al J 120061) . 
The imaging camera is described in lGunn et alj il998|). The pho- 
tometric system and calibration are described in 'Fukugita et aj 
( 1996), Hogg et al. (2001), Smith et al. (2002), Ivezic et al. (2004 
and T ucker et al.. ( 2006). The astrom etric calibratio n is in Pier et aU 
(2 003} ) and the data pipelines are in iLupton et alj EoOl) . LuptonI 
( 200a) for the deblender, Frieman et al. and Schlegel et al. (in 
preparation). 

The targeting pipeline cho oses targets for spectr oscopy from 
the imaging. A tiling algorithm iBIanton et all200'33) then assigns 
the spectroscopic fibres to the targets, the main source of incom- 
pleteness being the minimum distance of 55 arcsec between the 
fibres. This causes about 6% of galaxies to be missed, and these 
will be biased towards regions with a high surface density of galax- 
ies. The algorithm gives a more uniform completeness on the sky 
than a uniform tiling by taking into account large scale structure, 
but some bias is still present. 

The SDSS galaxies with spectra consist of a 'Main', flux- 
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limite d sample, with a median redshift of 0.104 iStrauss et alJ 
l2002h . a luminous red gala xy sample (LRG), approximately 
volume-limited to z ~ 0.4 ( Eisenstein et al. 200?) and a quasar 
sample ( Richards et al. 2002). The limiting magnitude for the Main 
spectra is r < 17.77, which is substantially brighter than that for 
the imaging so the redshift completeness is almost 100%. A typical 
signal-to-noise value is > 4 per pixel and the spectral resolution is 
1800. The redshifts have an RMS accuracy of ±30 km s~^. 

This paper uses galaxies from the DR4 version of the 
VAGC, which is desc ribed at http://wassup.ph vsics.nvu.edu/vagc/ 
iBlanton et all <2005j) describes the DR2 (fAbazaiian et alJl2004 
version. The catalogue is a publicly- available set of FITS files con- 
taining the Princeton reductions (http://photo.astro.princeton.edu 
^http://spectro. astro.princeton.edu I of the raw SDSS data with nu- 
merous additional derived quantities and matches to other surveys. 
The Princeton reductions are designed to improve on the origi- 
nal SDSS pipeline reductions in the publicly available DR4 Cat- 
alogue Archive Server (CAS; http://cas.sdss.org/astro/en ( and are 
used throughout unless otherwise stated. The set of objects con- 
tained is designed to match the SDSS DR4 and contains galaxies 
which either match a slightly more inclusive version of the Main 
sample criteria, are within 2 arcsec of a Main, LRG or quasar tar- 
get from the version of the photometry used for the targeting, or are 
within 2 arcsec of a hole drilled in an SDSS spectroscopic plate. 
The catalogue provides 6,851 square degrees of imaging coverage 
and 4,681 square degrees of spectroscopic coverage. The raw cata- 
logue contains 1,223,536 objects, 722,866 of which have spectra. 

The VAGC units differ from the CAS data in some quantities: 
the object brightnesses are given as fluxes in nanomaggies (where 
1 maggie = mag arcsec"^); the errors are inverse variance and 
the Petrosian radii are in pixels. These are converted to the units 
used in which the Main Galaxy Sample is defined, i.e. magnitudes, 
la errors and arcsec. The pixels are converted to arcsec using the 
pixscale values given in the full VAGC calibObj photometry out- 
puts. These are all close to 0.396 arcsec pixeP^. 

From the catalogue, the main file containing the imaging data, 
object_sdss_imaging.fits, was chosen as a base, and the rest of the 
required parameters for each object were obtained by matching the 
other datafiles to this one using a method appropriate for each pa- 
rameter and file in question. The matches are described further be- 
low. 

The magnitudes are corrected for galactic redden- 
ing using the usual correc ti ons de rived from the maps of 
ISchlegel. Finkbeiner & Davi3 il998l) . These are of order 
—0.1 mag, with a mean value of —0.14 mag and la vari- 
ation of 0.11 mag. The minimum and maximum values are 
—0.02 mag and —1.04 mag respectively, both in the u band. 

K-corre ctions are app l ied usi ng version 4.1.14 of the code 
described in Blanton et a l. 1 2003b). This provides a separate K- 
correction for each individual galaxy, but does not take into account 
galaxy evolution or dust. The K-corrections were made to a band- 
shift corresponding to a redshift of 0.1. Most of the datasets have 
approximately this mean. One could attempt to K-correct to the 
band-shift for the mean of each dataset, however, this would affect 
the dataset chosen and so this is not attempted here. 

2.1 Galaxy Samples and Properties 

The SDSS DR4 imaging outputs were used directly as supplied in 
the VAGC file object_sdss_imaging.fits and cross-matched from the 
full outputs available in the calibObj files. The DR4 spectroscopic 
outputs were used similarly from object_sdss_spectro.fits. 



The absolute magnitude range used for each bivariate LF was 
—24 < Mr — 5\ogh < —15 in 18 bins of 0.5 mag, with two 
exceptions: the JPG morphological type (see below) was binned 
in 12 bins of 0.5 mag for -24 < Mr - 51og/i < -18, and the 
environmental density was in 22 bins for —24.03 < Mr — 51og/i < 
— 19.63, due to its restricted redshift range (see below). Details of 
each bivariate LF are given in TableQ 

The upper limit of approximately 40 evenly spaced bins in 
the second parameter was chosen to give a reasonable execution 
time. The extent of the bins is based on visual inspection of the 
distributions of the input parameters, but the value of the LF in 
each bin is independent of the others. 

Our analysis allows for the sampling rate of the dataset, which 
is the ratio of the number of galaxies which have spectra taken to 
those meeting the criteria for spectra to be taken. The complete- 
ness maps as a function of imaging-only parameter pairs are given 
in Sj4| The overall completeness is estimated to be 93%. 6% of 
galaxies are missed due to the minimum 55 arcsec separation be- 
tween adjacent spectroscopic fibres on a plate, and the Main Galaxy 
Sample is about 99% complete overall ( Strauss et al. 2002). How- 
ever, this completeness is in terms of those objects that are tar- 
geted for spectroscopy in the first place, and in particular is sub- 
ject to an explicit Petrosian 50% light sur face brightness cut of 
/^app < 24.5 mag arcsec"^ for all objects. Blanton et al. ( 2005j) 
show that the faint end is incomplete. We therefore calculate the 
sampling rate for each bin in the bivariate LF and correct for this in 
the estimation of the LF (see fjs}. 

Because we calculate the sampling rate, the galaxies used in 
the process of obtaining our results do not necessarily have spectra. 
We therefore apply the set of recommended photometric flags de- 
scribed at http://www.sdss.org/dr4 to generate a clean sample. We 
also require the resolve_status flag to be survey .primary and the 
vagc.select flag to be main. 

For spectra, as mentioned, the VAGC Main-like sample cuts 
are slightly more inclusive than the SDSS Main sample, in partic- 
ular allowing in some binary stars by the removal of the cut for 
small bright objects . Here we reapply the Main criteria as given in 
IStrauss et alJ J2002h . For objects with spectra we also require the 
specprimary flag to be set, the primtarget flag to be set to either 
galaxy, galaxy _big, or galaxy _bright_core, the progname flag to be 
main, the zwarning flag to be and the platequality flag to be good. 
The VAGC does not contain a redshift confidence parameter equiv- 
alent to the zConf in the CAS, but bad spectra and redshifts are 
excluded by zwarning. 

In the optical, many spiral galaxies that are edge-on or near to 
this suffer internal extinction and reddening due to dust. Here we 
attempt to minimize this effect by restricting the sample to galax- 
ies with an axis ratio of less than that for an E7 elliptical galaxy. 
Galaxies with larger axis ratios than this are likely to be edge-on or 
close to edge-on spirals. Because the radial light profiles available 
from the SDSS and VAGC are either axisymmetric or specific to 
certain types of galaxies (de Vaucouleurs or exponential), we use 
the 25 mag arcsec^^ isophotal major and minor axes isoA and 
isoB in r and restrict the sample so that isoA/i soB < 10/3. This 
exclud es 30,407 (2.5%) of the galaxies. Fig. 4 of lVincent & RvdenI 
i2005h shows that the vast majority of the bright galaxies removed 
are those with an exponential or mixed exponential and de Vau- 
couleurs profile. Very few galaxies with a pure de Vaucouleurs pro- 
file are removed. At faint magnitudes similar numbers of galax- 
ies are removed for both profile types, so many more faint ellipti- 
cals are removed than bright. Here 'bright' depends on the profile 
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type, varying from A/,. — 51og/i < —21.1 for de Vaucouleurs to 
Mr — 51og/i < —19.7 for exponential. 

A potential alternative measure is the axis ratio derived 
from the adaptive moments I Stoughton et al. 2002). However, 
iKuehn & R vden l20o3) show that this would be affected by seeing 
for most of the galaxies in our sample due to their small sizes and 
consequent poor resolution. As our purpose is simply to exclude 
very elongated objects we use the isophotal cut. A similar example 
is lZibetti e t al. ( 200^, who require isoB/isoA < 0.25 in g, r and 
i, and a > 10 arcsec for their sample of edge-on galaxies. 

The VAGC also contains eyeball quality checks of spectra for 
17,422 objects, mostly bright (Mr — 51og/i < —15) or of low 
redshift (z < 0.01). The sample is not complete in any particular 
sense, but it does exclude numerous spurious objects such as bad 
deblends of large galaxies, HII regions, clearly incorrect classifi- 
cations and so on. We require that the quality flag be set to done 
with no other bit set, unless use.anyway is set. The requirement 
noticeably reduces the numbers of objects in some of the sparsely 
populated outlying bins of the LF. 

Besides the VAGC, we use data from five other sources: 
1) the Japan Participation Group (JPG) catalogue of eyeball- 
classified gal axy morphologie s, based on the SDSS Early Data Re- 
lease (EDR, IStoughto n et alj| 2002l). 2) the neural network mor- 
phological tvpes of lBalletalJ l20o4\ . 3) the public SDSS DR4 
CAS, 4) the Max Planck/Johns Hopkins (MPA/JHU) catalogues 
for galaxies and AGN at http://www.mpa-garching.mpg.de/SDSS 
and 5) the Pittsburgh-Carnegie Mellon Value-Added Catalogue 
(VAC; http://nvogre.phyast.pitt.edu/dr_value_added ( for galaxies. 
The CAS gives the eClass spectral type, the MPA/JHU the stellar 
mass and the VAC the environmental density. 

2.2 Petrosian and Sersic Magnitudes 

Petrosian magnitudes measure a constant fraction of the total light, 
in a model-independent manner. They are available i n the SDSS 
in a modified form from that introduced bv Petrosiar] il976h . The 
Petrosian flux is given by 

(■JVprp 

Fp = 2nr'dr'l{r') (3) 

Jo 

where rp is the Petrosian radius, which is the value at which the 
Petrosian ratio of surface brightnesses 

/o'gf27rr'dr7(r')/[7r(1.25' -0.8^)] 
^ J^2nr'dr'I{r')/{7vr^) ^"^^ 

has a certain value, chosen in the SDSS to be 0.2. The number A^p 
of Petrosian radii within which the flux i s measured is equal to 2 
in th e SDSS. Furth er details are g iven in lLupton. Gunn & Szalavl 
il999.) and Stoughton et alJj2002K . 

The Petrosian radii are not seeing corrected, which causes the 
surface brightness and concentration to be underestimated for ob- 
jects of size comparable to the PSF. However the seeing effect is not 
yet quantified, and there are other approximations such as using a 
circular as opposed to elliptic al aperture and not correcting for dust 
obscuration. As described in lstoughton e t al. (2002) the Petrosian 
aperture is not missing much flux compared to an ideal galaxy light 
profile, the amount missing being about 20% for a de Vaucouleurs 
profile and only 1% for an exponential profile. The effect of seeing, 
which would make the profiles tend towards a PSF, for which 5% 
of the light is lost, is also small for galaxies in the main sample. 
iBlanton et alj i2003d) show that the resulting luminosity density in 



r is very similar to that from the Sersic profile, the difference be- 
ing jo. 1 (Sersic) = jo. i (Petrosian) — 0.03. They suggest that the 
similarity shows that the true luminosity density is also of a similar 
value. 

The VAGC also contains Ser s ic fits to the galaxy light p rofiles. 
The Sersic profile ( Sersidll968t iGraham & Driverll2 00j) is ob- 
tained by generalising the de Vaucouleurs profile Tde Vaucouleurl 
1 19481) to have index i instead of i , giving 

/(r) = Joexp{-6„[(r/r,)i/"]}, (5) 

where b„ is such that half the total luminosity is wi thin Ve- The 
fitting procedure used in the VAGC is described in Blanton et al 
(2005a) and in more detail in the Appendix of Blanton et al 
( 2005b|). The axisymmetric Sersic profile is fitted to the azimuthally 
averaged galaxy light profile available from the SDSS database. 
The profiles are corrected for seeing, which is modelled using three 
axisymmetric Gaussians. The de Vaucouleurs and exponential pro- 
files corre spond to w = 4 and n = 1 respectively. 

Grah am et alj i2005h show the discrepancy between the Pet- 
rosian and Sersic magnitudes. In particular, the difference depends 
on the shape of the galaxy light profile. Their table 1 for the 
SDSS Petrosian aperture shows that for inverse concentration index 
(C/inv; ^2.4> . the difference is negligible above C/inv = 0.35, but 
below this increases to, where is the total effective half light sur- 
face brightness, fie — Mapp = 0.52, 1.09 and 1.74 mag arcsec"^ 
for C/inv ~ 0.30, 0.28 and 0.27 respectively for mean half- 
light surface brightness, and mpctrosian — Jnscrsic ~ 0.22 mag, 
0.38 mag and 0.54 mag in magnitude over the same range. Thus 
only a few percent of our galaxies are affected significantly. Their 
fig. 6 shows that the VAGC Sersic fits are consistent with the ex- 
pected differences. The C/inv of 0.28 corresponds to a Sersic index 
of 6. 

Although the Sersic index is an improvement in that it fits 
a more accurate radial light profile to a galaxy and is seeing- 
corrected, there are still some biases present. In their fig. 9, 
iBlanton et alJ i2005A show the residuals of their Sersic fits for 
1,200 simulated galaxies. As a function of n, the value of Uout — 
riin, where Uout is the fitted value of n, decreases monotonically 
from less than -0.05 at n = 1 to -0.1, -0.25, -0.5 and -0.65 at 
n = 2,3,4 and 4.5 respectively. Thus a de Vaucouleurs galaxy 
(n = 4) is assigned an index of n ~ 3.5 and the indices are all 
systematically underestimated. The error of 0.5 is comparable to 
the bin size in the same area of our bivariate luminosity-Sersic in- 
dex distribution. The range of the error as shown by the quartiles 
also broadens from around 0. 1 to 0.5 over the same range. The fit- 
ted 50% light radius and flux decrease similarly over their range to 
around 0.8 and 0.9 of their true values. Also, the fitted radius and 
flux decrease with n, the fitted flux with radius and the fitted ra- 
dius with flux. However, the fitted n does not decrease with radius 
and flux. The overall performance for n is characterized as good 
because the bias is similar to the uncertainty and comparable but 
opposite to that from the assumption of axial symmetry. The over- 
all fits are thought to be approximately correct for Sersic-shaped 
galaxies and supply a seeing-corrected estimate of size and con- 
centration for the others. 

Here we investigated the LFs for both Petrosian and Sersic 
profiles. The LFs are generally similar and because the samples 
used are based on the SDSS Main Galaxy Sample, which is based 
on Petrosian magnitudes, we present the results based on those 
magnitudes. The derivation of a fully Sersic-based sample is be- 
yond the scope of this paper. 
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Table 1. Parameters and ranges used in the bivariate luminosity distributions. See text for further details on the parameters and data sources. All data are 
based on the SDSS DR4, with the exception of the JPG catalogue, which is from the EDR. Ng^i is the number of galaxies, to the flux limit, with spectra and 
selected for the luminosity function respectively. The overall sampling rate is Ng^i^p^^j^,^ /^galiinaging ■ ^'-^ ''^^ best-fitting faint-end slope, 

characteristic magnitude and normalisation respectively for the overall best-fit Schechter function for each LF. 
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2.3 Neural Network Morphological Types 

1,875 SDSS galaxies to r < 15.9 in the SDSS EDR have be en clas- 
sified into morphological types by Nakamura et al.' ("2003*), form- 
ing the JPG catalo gue. The system used was a modified version of 
the T-type system ide VaucouleurjlT95 91 with the types being as- 
signed in steps of 0.5. The corresponding Hubble types are E=0, 
S0=1, Sa=2, Sb=3, Sc=4, Sd=5 and Im=6. 

We have prev iously shown jBalLeLalJ2QQ4|) that artificial neu- 
ral networks (e.g. lBishodll995l : lLahav et alJll993) are able to as- 
sign types to galaxies in the SDSS with an RMS accuracy of 0.5 on 
this scale, using the JPG catalogue as a training set. This is the same 
as the spread between the human-classified types by the members 
of the JPG team. 

We used the same procedures here, updated for DR4 over 
DRl and with some modifications. The DR4 CAS was matched 
to the JPG catalogue using a tolerance of 2 arcsec and a flux limit 
of r < 15.9. Unassigned types (—1) and galaxies flagged as be- 
ing likely to have bad photometry were removed as before. The 
clean photometry flags described above were also applied here. Af- 
ter these cuts 1,290 galaxies remained, 1,131 with spectroscopic 
redshifts. The training and test sets were 645 each with no overlap, 
from which, as before, galaxies with severely outlying parameters 
and targets (> IOct from the mean value for the parameter, gen- 
erally indicative of a measurement error) were iteratively removed 
for each parameter in turn. The network trained on the 645 was 
applied to the DR4 sample using the CAS photometry. One could 
also retrain on all 1 ,290 for the final network to be applied to the 
sample, but it would make very little difference. 

It was also found th at when trained on the full set of 29 param- 
eters in 'Ball et alj i2004h the types were biased towards early type 
with increasing redshift to a greater extent than in the JPG cata- 
logue. However, when one trains on the purely morphological sub- 
set of the parameters (all those except the magnitudes and colours), 
the bias is similar to the JPG, reflecting the flux-limited nature of 
the sample, but the RMS variation between network type and target 
type is not significantly larger Hence the latter set is used as the 
training set here. 

The network also has a tendency to 'avoid the ends of the 
scale', resulting in few very late type galaxies and a slight bias 
away f rom early type galaxies. This is shown in Fig. 1 of lBall et alJ 
yOOJ). This is likely due to the severe lack of galaxies in our train- 
ing set of type T = 5 or later compared to the earlier types, and the 
importance of the concentration index in the morphological train- 



ing set, which shows a similar deviation away from the ends of the 
scale for both early and late types. 

We restrict the absolute magnitude range of the morphologies 
presented to Mr — 51og/i < —18. This is due to a lack of galaxies 
in the training data at fainter magnitudes than this. Thus the neural 
network types are not extrapolated from the training set in apparent 
or absolute magnitude. 



2.4 SDSS DR4 VAGC Data 

The bivariate LF parameters described in this section are the in- 
verse concentration index, Sersic index, absolute effective surface 
brightness, reference-frame colours and absolute radius. The mor- 
phological parameters are measured in the r band, since this band 
is used to define the aperture through which Petrosian flux is mea- 
sured for all five bands. 

The inverse concentration index C/inv is -Rso / Rf>o where iiso 
and i?go are the radii within which 50 and 90 per cent of the Pet- 
rosian flux is received. The inverse is used because it has the range 
0-1. 

The Sersic index n is monotonically related to the concen- 
tration (e.g. Graham et al. 2005), if the latter is measured using 
Sersic radii. Here we use Petrosian circular radii and for the bi- 
variate luminosity- Sersic index distribution we use log 71, as the 
index is in the exponent in the equation defining the Sersic profile 
(equation |5j. 

The absolute effective surface brightness used here is given by 

A^Rso =mr + 2.51og {2itRIo) - 101og(l + z)-K, (6) 

where m,. is the r band magnitude and K is the K-correction. The 
evolutionary correction term is set to zero. 

The reference frame colours are taken from the K-corrected 
values described above. All ten colours from u, g, r, i and z were 
investigated. 

The absolute radius of a galaxy at redshift Zgai is calculated 
from the apparent Petrosian radius using the f^mattcr — 0.3, — 
0.7, h = 1 cosmology. It is also shown logarithmically and is given 
by 

^90 d{z) 

-Kabs — — - — [ , (') 

1 + Z 
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where d is the relativistic coordinate distance 



^0 Jo 



-gal 



(8) 



2.5 Outputs from additional SDSS Catalogues 

The additional catalogues used are all based on SDSS DR4 and use 
the same Euclidean space, fimattcr = 0.3, JIa = 0.7, as used here. 
The quantities all require spectra to be calculated and are therefore 
matched using the unique spectroscopic identification, both in ob- 
ject and time of observation, provided by the quantities plate, mjd 

and fiberid. 

[ ■ eClass jConnoUv et 111119951: Iconnollv & Szalavlll999l: 
lYip et alj206i) is a continuous one-parameter type assigned from 
the projection of the first three principal components (PCs) of the 
ensemble of SDSS galaxy spectra. The locus of points forms an ap- 
proximately one dimensional curve in the volume of PC 1 , PC2 and 
PC3. This is a generalization of the mixing angle (j) in PCI and PC2 



tan 



(9) 



where ai and a2 are the eigencoefficients of PCI and PC2. The 
range is from approximately —0.6, corresponding to early type 
galaxies, to 1, late type. The eClass is also robust to missing data 
in the spectra used for its derivation, and is almost independent of 
redshift. The quantity is not given in the VAGC and so was matched 
to the public SDSS DR4 data from the CAS. 

The MPA/JHU value-added catalogues at 
Phttp://www.mpa-garching.mpg.de/SDSS include the quantity 
rml50, the logarithm of the median dust-corrected stellar mass 
in solar units. It is calculated from the stellar mass to light ratio 
predicted by a large li brary of models of star f ormation history and 
is described further in lKauffmaim et ai]i2003h . 

The VAC contains a measure of galaxy environment as the 
parameter density_z005095. This is the distance to the A''th near- 
est neighbour within ±1000 km in redshift for 0.053 < z < 
0.093. The ±1000 km is used to minimize contamination from 
interlopers. Galaxies for which the survey edge is reached before 
the A'^th nearest neighbour are excluded to avoid downward bias 
in the estimated densities near the survey edge s. Here the value 
of A'^ used is 5, following [Balog h et alJ i2004) who choose this 
value to approximate iDressleJ ^1980^ who uses A = 10 before 
correction for superimposed galaxies. The surface density of galax- 
ies is then given by Ejv = N/Tvd%. The nearest neighbour must 
be of magnitude r < 17.7 at z — 0.093, which corresponds 
to Mr — Slogft < —19.63. This cut ensures a uniform density 
measurement over the redshift range of the sample rather than the 
density decreasing with redshift due to missing faint galaxies. The 
value of h used in preparing this catalogue was h = 0.7 so the 
densities are given here in units of Mpc"^. 



3 BIVARIATE LUMINOSITY FUNCTIONS 

Unlike the monovariate LF in which the Schechter function is of- 
ten used, there is no common functional form for the bivariate 
LF, so a nonparametric estimator is chosen. We use the popular 
stepwise maximum likelihood (SWML) method of lEfstathiou et alj 



stepwi ! 
h988h 



following the extension to a bivariate distribution by 



defined error properties teendall & Stuartlll96llj and in the com- 
parison between LF estimators o f lTakeuchietalJj2000l) the SWML 
was shown to be a good estimator. 

The SWML gains independence of inhomogeneities in the 
galaxy distribution by assuming that this is the case via the univer- 
sal form: n{L, x) = (/)(L)p(x). This means that the ^(L) shape is 
independent of its normalisation, which then has to be found sepa- 
rately. The data must also be binned over the ranges of the param- 
eters considered. 

Here the bivariate LF is given by tp{L, X) or 4>{M, X) where 
X is the second parameter. 

For a galaxy which is observable if its luminosity L lies in the 
range I/min-^max and X lies in the range Amin-Amax (these lim- 
its in general being redshift dependent), the probability of seeing 
that galaxy with luminosity Li and A = A^ at redshift Zi is 



<x^{U,X,) f{L„X,) 



xjj{L,X) f(L,X) dLdX, (10) 



where / is the completeness function, calculated for each bin in 
the bivariate LF. This automatically takes into account both the 
6% sampling incompleteness resulting from galaxies more closely 
spaced than the spectroscopic fibres, and the further 1% from ob- 
jects missed in the Main Galaxy Sample (see i|2.H . It does not take 
into account any incompleteness in the imaging itself, which is dis- 
cussed further in Appendix]^ below. The A limits are a function 
of z if the sample is explicitly or implicitly selected on A. The 
samples are not explicitly selected on A if the quantity requires a 
spectrum. 

The likelihood C — 11 P» then maximised with respect to 
tp{L,X). In practice, one maximises the log-likelihood In £ = 
Eln Pi- 

For the SWML method, the discrete version of this is used. 
ip{M, A) is parametrized as the number density of galaxies 

i/'(M,A)=V,fe {j ^1...Nm; k = l...Nx) (11) 

in Nm and Nx evenly spaced bins in absolute magnitude M~ < 

M, < M+ and A^ < A^ < A+, where 



and 



Mt = Mi ± 



Ai. — Afc ± 



AM 



AA 



(12) 



(13) 



ISodre & Lahavl Jl993h . The maximum likelihood method has well 



Magnitude bins are preferred to luminosity bins as the distribution 
is more evenly spread. 

The likelihood is then given using the discrete equivalent of 
eguation llOl 

JVg Nm Nx 

In ^ = E E E ^'^^ In t^^-fc fiM.,Xi)] 

Ng /Nm Nx \ 

— E In I E E ) + constant, (14) 

i=i \j=i fe=i / 

where 

C 1 if - AM/2 ^ Mi - Mj < AM/2 
W,jk = 1 and - AA/2 A, - A^ < AA/2, (15) 
I otherwise 



© 2006 RAS, MNRAS QOQ.mEU 



Bivariate Galaxy Luminosity Functions in the SDSS 7 



and 



Hijk — 



1 



AMAX 



M' pX' 



f{M,X) dM dX, (16) 



where A/' = max[M" , min(M+, Mi^„J] and X' = 
max[X- , min(X+, X/i„J]. 

The constant in <14> is fixed by using a Lagrangian muhipher 
A: the constraint 



g = J2Yl ipjkAMAX -1=0 

3 k 



(17) 



is applied and the new likelihood ln£ = ln£, + Xg^ijjjk) is max- 
imised with respect to ipjk and A, requiring A = 0. 

The maximum likelihood estimate d\nC /dtpjk = is then 



Ellfi w,,k 



(18) 



where ipim is from the previous iteration. 

Following Efstathiou et al., the errors on the parameters are 
estimated using the fact that the MLM estimates ipjk are asymptot- 
ically normally distributed with the covariance matrix 



COv(V'jfc) = I ^(^jfe), 



(19) 



where I is the information matrix (their equation 2.13b). 

For the normalisation, the estimator of the space density n of 
galaxies is 



/n = ^ w{zi 



S{z)w{z) dV. (20) 



The selection function S from Li to L2 is 

^2 "^2 / -^2 -^2 

S{z)^Y.Y.^{L,X) Y.Y.^{L,X), (21) 

where Li = max(Lmin(2), Li), L2 = min(Lmax(2:), L2) and 
similarly for X. The limits on the sums in the numerator depend on 
the K-correction, which is determined using the average SED for 
the sample. These limits will in general include partial bins with 
appropriate weighting due to the summation limits not falling at 
the edge of a bin. 

We use the weighting function 



1 



\l + A-KfnJs{r,)S{z)y 



where 



Mr.) = 



^i{r)dr, 



(22) 



(23) 



r is the distance and ^(r) the real space galaxy two point correlation 
function. 

The variance in the correlation function can be estimated us- 
ing a complex formula which minimizes the variance in the func- 
tion and involves the three and four point correlation functions. 
These are more difficult to measure that the two point and require a 
large sample such as the SDSS. The weight function in eQuation l22l 
was found to be a good approximation by Hamilton ( 1993) and is 
widely used. The function minimizes the variance in the estimate of 
the number density so long as is much less than the survey depth. 
Here 47rJ3 = 32000 Mpc^, using the ^(r) from Zehavi et al. 
j2002h which is eiven as f fr) = (r/6.1±0.2 ft"^ Mpc)"^ '^^±" "^ 



for0.1/i~^ Mpc ^ r ^ lQh~^ Mpc. The effect of using an incor- 
rect value would be to change the nor malisation but not the shape 
of the LF. For example. lLovedavl fe004^ uses the same values of J3 
and ^ as here and finds that halving J3 reduces the estimated den - 
sity by 7%. An update to lZehavi et £d]j2002h . lZehavi et aljf2004bh . 
has a similar ^ value to that used here. 

Equation 1201 is solved iteratively as it contains w, which 
contains n. The minimum-variance w(z) is described further in 
iDavis & Huchrai 1.1982.) and .Hamiltoa jl993h . 



4 RESULTS 

4.1 Morphological Type 

The LF bivariate with the JPG catalogue morphological type T!ipg 
obtained here is shown in Fig. Q and that with the ANN morpho- 
logical type Tann in Fig.|2| The corresponding completeness maps 
(number of galaxies with spectra / number of galaxies, in each bin) 
are shown in Figs. lAl| and lA2l There are no obvious biases within 
the regime for which galaxies are present, the completeness being 
consistent with 90% for the full range of apparent magnitude for 
TjpG and 60% for Tann, the latter being expected from the use of 
the whole survey area which is only about 70% complete in spec- 
troscopy (Appendix!^. Fig. lA2l again shows the lack of types later 
than 4.5 assigned by the ANN (see i|2.3> . If the training set with- 
out the clean photometry flags were used, hence leaving in more 
late type galaxies, the ^j) greyscale would extend further right by 
approximately the width of one of the bins, so the LF for very late 
types would still not seen. 

The bivariate luminosity-morphology distribution was com- 
puted for the same JP G catalogue and classification scheme by 
iNakamura et alj i2003h . However, whilst being drawn from the 
same catalogue there are some differences in the method, leading 
to an altered final sample for the LF. A different subset of the EDR 
galaxies is used here, with 1,150 instead of 1,482; the sample is 
divided into smaller bins and t he normalisation is without the ex- 
plicit redshift cut used in Naka mura et alji2003h . The smaller value 
is 1,150 due to the masking applied to the datasets here. Of these, 
1,019 have redshifts and 998 are selected for the LF. 

Naka mura et alJ show that the overall LF is well fit by a 
Schechter function with parameters a = —1.10, M* = —20.65 
and (f)* = 0.0143. We obtain a = -1.30, M* = -20.58 and 
0* = 0.0098. Their LF is divided into the ranges < Tjpg < 1.5, 
1.5 < Tjpg < 3, 3.5 < Tjpg < 5 and 5.5 < Tjpg < 6. 
For 0-1 and 3.5-5 the faint-end slope is shallower, at a = —0.83 
and a = —0.71. 1.5-3 is steeper, at a = —1.15. The early types 
are also brighter, as expected, being M* — —20.75 compared to 
Al* = —20.30 for types later than 1. In our LFs there is no clear 
trend between the types, although the binning is narrower and the 
earliest type bin appears brighter and with a shallower faint-end 
slope as expected. 

In the LF split by Mr, the main trend is for higher normali- 
sation at fainter magnitudes, as expected for increased numbers of 
fainter galaxies. For bright galaxies, the normalisation with Tjpg 
increases from (f> ~ 10"'* to ~ 5 x 10~^ between types 1 and 
but otherwise is consistent with a flat distribution across the range 
in Tjpg at all magnitudes, with a possible decrease at types later 
than Tjpg = 4. 

The neural network types enable the LF sample to be en- 
larged from 998 to 20,891 galaxies. This is the first time the LF 
with morphology assigned in this way to the resolution of indi- 



© 2006 RAS, MNRAS 000.1711231 



8 A^. M. Ball et al. 



vidual Hubble types has been calculated for a sample of galax- 
ies of this size. The trends hinted at in Fig. Q become clearer, 
with the types earlier than Tann = 0.5 now showing a declin- 
ing faint-end slope. The 0.5 < Tann < 1-5 are as bright but more 
abundant at fainter magnitudes. The spirals from Tann = 1.5 to 
Tann = 3.5, corresponding to SOa to Sbc, then show no differ- 
ence in their LFs. The late types to 4.5 are then slightly fainter 
again, with a slightly steeper faint-end slope. Each slice appears 
consistent with a Schechter function. For the LF split by Mr, the 
galaxies brighter than Mr — 51og/i = —21 are now clearly more 
abundant at Tann < 1, with a weaker trend seen for —21 < 
Mr - 51og/i < -20. The rest to Mr - 51og/i = -18 are still con- 
sistent with a flat distribution, except at Tann < 0.5 where they 
decline. The changes at the early type end are at least as strong as 
those indicated, as the ANN bias here is, as mentioned, away from 
early types and not towards. There are very few galaxies later than 
Tann = 2.5 and brighter than Mr — 51ogft = —22. 

In |de LapparenH |2003), it is suggested that cross- 
contamination between Hubble types will make the LFs tend 
towards Schechter functions even if the underlying LFs for the 
giant (as opposed to dwarf) galaxies are Gaussian. Here the 
—0.5 < Tann < 0.5 bin in particular could be suffering from 
this effect if the true distribution were Gaussian brighter than 
Mr — 51ogft = —20 a nd the fai nter magnitudes were contami- 
nated by dwarfs. Nakam ura et alj (2003) do not find this to be a 
significant effect in their results as their fainter galaxies do not 
have the softer cores of late r type dwarfs. Howev er, the Gaussian 
LF is supported by fig. 1 of lBemardi et alJ j2003), who show LFs 
for ellipticals with a turnover at Mr — 51og/i ~ —20.5. The LFs 
for the later types will similarly contain dwarf and giant galaxies 
and will therefore also appear Schechter-like. 

Several recent papers discuss blue dwarf spheroidal galax- 
ies. These have the same concentrations and smooth appearance 
as dwarf ellipticals, but are much bluer in colour. Whilst much 
more prevalent at higher redshifts (e.g. Ilbert et al. 2006), for lo- 
cal galaxies these objects are discussed by Driver et al. (2006). In 
their fig. 6, they show that when the blue spheroids are removed, 
the E/SO LF gains a much more rapidly declining faint-end slope, 
with Q changing from -0.9 to -0.4. This could be addressed in our 
work by adding colour as a classification criterion, moving towards 
a spectro-morphological classification. 

It may also be the case that the morphological types are not 
well defined below a resolution of about 0.5 T types and will there- 
fore inevitably be cross-contaminated. However with the large sam- 
ple here the LF differences are significant, as indicated by the sizes 
of the error bars. 

Another recent determination of t he LF bivariate with Hubble 
type is presented in Read & Trenthan] i2005h . This shows LFs fol- 
lowing the same trends as seen here, and with only those for types 
Sa and Sb being similar. 



4.2 Inverse Concentration Index 

Various papers (e.g. lShimasaku et all200lh have shown that the in- 
verse concentration index R50/R90 (C/inv) is correlated with the 
morphological type, and indeed can be used as a simple classifier. 
Here we show the LF bivariate with concentration index in Pet- 
rosi an radii, using the same R50 / Rgo ratio. 

IShimasaku et alj show that, whereas the C/inv does correlate 
with morphology, the correlation is not perfect. In particular, for 
types E and SO, the index varies from around 0.3 to 0.35, whereas 
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Figure 1. LF bivariate with JPG catalogue morphological type TjpG for 
998 galaxies. The contours show the value of logjQi/)(M, X), the bivari- 
ate luminosity-X distribution in each bin. The second panel shows vertical 
slices of ^p{M,X) subdivided by X for the ranges given and the lower 
panel shows the horizontal slices '4>{M, X) subdivided by Mr- Although 
the contours show all galaxies, a point appears for the slice only if created 
by five or more galaxies, as below this threshold Poisson fluctuations dom- 
inate. In the rest of the LF plots (Figs. l2l[T3l this is increased to twenty (see 
text). Here X is Tjpq. The errors are from the inverse of the information 
matrix (see text). Note that for clarity the axes in this figure and that for 
Tann differ from subsequent figures due to the brighter cutoff in absolute 
magnitude. 
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Figure 2. As Fig. but for the neural network morphological type. The 
trends in Fig.0become clearer due to the sample size of 20,891 compared 
to 998. 



to the prevalence of dwarf ellipticals in clusters and is discussed 
further below. 

Subdivided by Mr, at the bright end, -24 < Mr - Slogft < 
—21, the C/inv shows a peak at C/inv ~ 0.31. A pure de Vau- 
couleurs profile corresponds to C/inv = 0.3, therefore this peak 
corresponds to the population of ellipticals. The slight offset above 
0.31 may be due to the bias towards the C/iiw for a PSF, which 
is 0.5 (Gaussian), due to many of the galaxies in the sample being 
poorly resolved, in the sense that their isophotal axis ratio rather 
than their shape measure is used ( i|2. 1> . However, this region is also 
populated by SO galaxies. At -21 < Mr - 51og/i < -19.5, the 
peak at low C/inv becomes broader and less pronounced, increas- 
ing to C/inv ~ 0.33. At -19.5 < Mr - 51og/i < -18 the distri- 
bution becomes more symmetrical, peaking around C/inv ~ 0.38. 
At —18 < Mr — 51og/i < —16.5 the distribution is approximately 
symmetrical around C/inv = 0.43, which corresponds to an ex- 
ponential profile. Besides the bright peak, the distributions are not 
otherwise obviously bimodal, just broad. 

The completeness increases from ~ 40% to 50% over the 
range 0.3 < C/inv < 0.4 and remains at this level until C/inv ~ 
0.55. This may reflect the bias due to the closest possible spacing of 
55 arcsec for SDSS spectroscopic fibres (32)' more galaxies in en- 
vironmentally dense regions, which tend to be early type and there- 
fore of C/inv in this range, will be missed. If this is the case then 
the peak in low C/inv seen should be slightly more pronounced 
than shown. 

iBlanton et alj i200lh present the LF bivariate with the Pet- 
rosian C/inv for the SDSS EDR. The same broads trends as seen 
here are present: a peak near C/inv = 0.3 corresponding to galax- 
ies with a de Vaucouleurs profile, particularly evident for the bright- 
est galaxies at —23.5 < Mr — 51og/i < —22, then becom- 
ing less prominent as the distribution broadens to spread as far 
as C/inv ~ 0.5 at 20.5 < Mr - 51og/i < -19.0 and becom- 
ing symmetrical about Clinv ~ 0.43, the exponential profile, at 
Mr — 51og/i < —17.5. The de Vaucouleurs peak is offset slightly 
above 0.3 in a similar way to our results. Unlike our LF, the nor- 
malisation split by Mr increases then decreases, whereas ours in- 
creases monotonically from bright to faint, presumably due to a 
different method of normalisation. Subdivided by C/inv, the LFs 
show the usual trend of bright and shallow faint-end to dim and 
steep faint-end. Compared to our results which cut in axis ratio, 
there is no obvious difference in the widths of the distributions in 

C/i nv ■ 

iNakamura et alJ i2003h also present the bivariate luminosity- 
C/inv distribution for their JPG sample described in i|4. II dividing 
the sample at C/inv ~ 0.35 into early and late types. As with our 
results above C/inv = 0.3, the LF does not show a marked decline 
at the faint end. They show that their sample is not dominated at 
the faint end (Mr }^ —19) by dwarf ellipticals, which have softer 



it increases to almost 0.5 to type Sc. It then stays in this range for 
even later types. 

The bivariate luminosity-C/inv distribution is shown in Fig. 
|3|and the completeness map in Fig. lA3l The expected broad trend 
of fainter galaxies being less concentrated is seen. Subdivided by 
C/inv the LF resembles a Schechter function with a steeper faint- 
end slope for less concentrated galaxies. For the bins outside 0.3 < 
C/inv < 0.6 the normalisation is lowered due to lack of galaxies. 
The completeness similarly becomes very low and noisy outside 
this range. For 0.3 < C/inv < 0.4 the LF shows an upturn at 
magnitudes fainter than Air — 51og/i = —18. This may be related 



4.3 Sersic Index 

The seeing-corrected Sersic index is mo notonically related t o the 
Sersic-radius-based concentration index iGraham et alJl20()^ and 
as such should show similar trends to those for the Petrosian-based 
C/inv presented in i)4.2l However, unlike our concentration index, 
the Sersic index is seeing-corrected. We plot the index logarithmi- 
cally as explained in i|2.4l The bivariate LF is shown in Fig.|4]and 
the completeness map in Fig. IA4I The indices again correspond 
to profile types, with log n = (n = 1) being exponential and 
log n = 0.6 (n = 4) de Vaucouleurs. The completeness decreases 
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Figure 3. LF bivariate with /Jso / i?90 inverse concentration index. The ex- 
pected trend of fainter galaxies being less concentrated is seen, but since 
the result is based on 251,744 galaxies (see TablefTl. much further detail is 
visible (see text). 



from around 55% at log n ~ —0.2 to 40% at log n = 0.65. Again 
this could be due to bias against early types due to fibre spacing. 

Subdivided by Mr, the distribution is clearly bimodal to 
fainter magnitudes than seen for C/inv, down to Mr — 51og/i = 
— 18. The high-n peak is around log n — 0.65 for —24 < 
Mr — 51ogft < —21, higher than de Vaucouleurs galaxies. For 
—21 < Mr — 51og/i < —19.5 the peak shifts lower, to log n ~ 
0.575. Below Mr — 51og/i — —19.5, the low-n peak becomes 
prominent and is centred at log n ~ 0.15. Again this is higher than 
the exponential profile of log n = 0. At this value, the number of 



galaxies has dropped sharply, although the completeness remains 
flat until log n ~ -0.2. 

Subdivided by n, the LFs are again Schechter-like and go from 
bright and shallow faint-end to dim and steep faint-end. As with 
C/inv, the faint end of the LF for high n, whilst decreasing, does 
not resemble a Gaussian. Thus if contamination by dwarfs affects 
the faint-end LF for morphology and C/inv, it is also at work for 
n. The LF turnup at Air — 51og/i — 18 seen for C/inv is less 
prominent here, merely hinted at for 0.3 < log n < 0.8. 

As described in i)2.4l the indices are systematically offset, so 
that a true de Vaucouleurs galaxy with log n = 0.6 is here as- 
signed log n ~ 0.55. However, the value of log n for a PSF is -0.3 
(n — 0.5), so this effect may again be at work, although the PSF is 
convolved in the Sersic model fitting so this seems unlikely. 

The bivariate lum inosity-log n distribution is also presented 
by iDriver et ail j200d) for the MGC. There it is subdivided into 
n ^ 2 and n < 2 (log n = 0.3) and fitted with a Schechter 
function. The respective parameters change from a = —1.25 to 
a = -0.66, M* = -19.48 to M* = -19.35 and (p* = 0.0129 
to (f}* — 0.0087. The high-n LF again is Schechter-like. The less 
obvious distinction in their two distributions than those seen in their 
LFs for the MGC continuum measure and colour is ascribed at least 
in part to dwarf ellipticals, which spread across the log n = 0.3 
boundary. As with our results from i|4. II colours are needed, which 
they go on to consider, finding a clear bimodality in the bivariate 
distribution of core ?t — r colour and n. This and their spread in 
Sabc morphologies between the two peaks is argued as evidence 
for bulges and discs being the fundamental components and is dis- 
cussed further in the context of our results below. 



4.4 Surface Brightness 

The bivariate luminosity-surface brightness (SB) distribution has 
been investigated by several previous authors and is also known as 
the bivariate brightness distribution (BBD). It is important as it al- 
lows one to better quantify sele ction effects in a surv ey. Here we 
compa re our results to those of iBlanton et alj i200ll).] i 



IShen et al 



I2OO3I) and Blanton et^ llOOSd) from the SDSS and lOriver etal 



(2005) from the Millennium Galaxy Catalogue. Further recent 
studies at low red shift in the optical are .de Jong & L acev (200^), 
ICross et alj j200lh and ICross & Driveil i2002h . The former stud- 
ies a more local sample of spirals in more detail and the latter 
two contain slightly earlier results than the MGC, from the 2dF- 
GRS. The SDSS studies use similar data and the same passbands 
as our data. The MGC represents the current state-of-the-art in 
the optical by probing to a deeper limiting surface brightness of 
26 mag arcsec"'^ in the B band over smaller but still substantial 
(37.5 deg2) area of sky. The broad trends here and in other re- 
cent work are in agreement with earlier work (e.g. iFreemanI 1 970l 
IPhillipps & Disnev. .1986^ etc.), but are seen at higher signal to 
noise. 

In the SDSS, the Petrosian and Sersic radii are derived from 
axisy mmetric profiles. As neither the Petrosian f i|2.4l Blanton et al] 
2003c) nor the Sersic apertures are missing much flux for most 
galaxies, the SBs will therefore be systematically overestimated as 
the semimajor axis of a half-light elliptical aperture is always at 
least as large. The MGC study uses elliptical apertures. 

The bivariate LF is shown in Fig.|5|and the completeness map 
for apparent surface brightness is shown in Fig. IA5I Also shown 
are contours at 10^-10^ Mpc"^ which show the regions in- 
side which these volumes were probed by the survey. Note that 
the axes of the top panel are extended relative to the LF slices 
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Figure 4. LF bivariate with seeing-corrected Sersic index n. The Gaussian- 
like LF for the galaxies of high, n, corresponding to early morphological 
type, is clearly seen. 



shown in the lower two panels. This shows that normalisation of 
LF drops before the volume becomes small due to limiting surface 
brightness and that the volume at Mr — 51og/i < —15, for which 
the LF is quoted, is at least of order 10^ h~'^ Mpc^. As a func- 
tion of /iapp, the completeness varies between 45% and 60% for 
19 < Mapp i$ 23 and outside this range becomes noisy and drops 
to being consistent with zero at /iapp ~ 17 and /iapp ~ 24. As a 
function of apparent magnitude, the completeness is between 55% 
and 60% in the range 12 < r < 16.5 and drops outside this to 
40% at r = 17.77, and low values at the bright end, consistent 
with zero brighter than r ~ 11.5. There is a clear lack of faint high 
surface brightness or bright low surface brightness galaxies. Fig. 



IA5l shows the incompleteness for spectra and, as described in Ap- 
pendix 1^ Blanton et al. ( 2005c) estimate the incompleteness for 
the photometry. Their ranges at which the completeness dropped 
below 90% were /iapp > 22, when both /iapp < 19 and r > 17, 
and /iapp ^ 18.7. They show that the SDSS galaxy distribution 
drops off well inside these completeness limits. Therefore, our dis- 
tribution will similarly and the overall LF is unaffected away from 
the bright and faint end. 

In the LF subdivided by Mr, we see the same overall distri- 
bution and trends as the previous studies. The /iRjg is a lognor- 
mal distribution with a roughly constant peak at bright magnitudes. 
This then broadens and dims at fainter magnitudes. The peaks are 
/iRso ^ 20.1, 20.5, 21.0, 21.5 and ~ 21.75 in steps of 1.5 magni- 
tudes for —22.5 < Mr — 51og/i < —15. No bimodality is evident, 
meaning that if the galaxies are divided into one or more popu- 
lations, as suggested by some of the other measures here such as 
colour, the surface brightness LF is insensitive to it, or the distri- 
butions overlap to such an extent that there is no central dip, with 
perhaps just a broadening. 

Subdivided by /irjq, the LFs appear Schechter-like, with a 
fainter M* and steeper faint-end slope towards fainter surface 
brightness, as expected, although the brightest bins lack galaxies. 
This is also consistent wi th previous results. 

iBlanton eta i'^2001') present the BBD for the SDSS EDR. Al- 
though the o veral l luminosi ty dens ity found by them has been su- 
perseded by iBlanton et alj i2003 ch . their BBD is consistent with 
ours. Petrosian magnitudes are used. They also use the EDR cal- 
ibrations of the passbands u* , g* , etc. but the difference here is 
small. They find for -23.5 < Mr* < -20.5 that the peak is at 
/iRjo ~ 20, the same as ours for —22.5 < Mr — 51og/i < —19.5. 
The distribution then broadens and the peak dims to fJ,R^g ~ 
20.5, 21.5 and 22 over the steps of 1.5 mag to Mr — 51og/i < 

— 16.0. Given the slightly different binning and the broadness of 
the peaks these are consistent with our values. At the bright end 
(Afr* < —19) the distribution drops off well above the EDR sur- 
face brightness cut of /^Rsq < 23.5 but at the faint end this may 
be an issue, as they are not able to s how th e upturn in the faint-end 
LF later reported bv IBlanton et alj j 2005^. We see similar results 
relative to the /tRgo < 24.5 cut. IBlanton et alJ also plot the LF 
versus SB, finding as we do that the LF appears Schechter-like, 
with a fainter M* and steeper faint-end slope towards fainter sur- 
face brightness. We see a clearer downturn in the LF in the range 
18.2 < MR.. <19.5. 

IShen et alJi2003h use data from a sample similar to the SDSS 
DRl. In their fig. 14, they show the surface brightness distributions 
for early and late type galaxies, defined by concentration index less 
than and greater than 2.86 respectively. The measures use Sersic 
magnitudes and h = 0.7, which we convert to h — 1 here. Late 
types show mean SBs of around 20.2 mag arcsec"^ for —21.2 < 
Mr — 51og/i < 20.2. This then drops to 20.5 mag arcsec"^ for 
Mr — 51og/i — —18.2 and more rapidly to 21.5 mag arcsec"^ for 
Mr — 51og/i = — 15.2, the faintest level. The early types are around 
20 mag arcsec"'^ for Mr — 51og/i — —21.2, 19.5 mag arcsec"'^ 
for Mr — 51og/i — —20.2 to 19 mag arcsec^^ for Mr — 51og/i — 

— 18.2, the faintest. The distributions again broaden at fainter mag- 
nitud es. These are approxi mately consistent with our results. 

IBlanton et alj <2005d) use the SDSS DR2 and show the over- 
all Petrosian SB distributions for -20.5 < Mr < -12.5. 
Over the range —20.5 < Air < —15 their peak varies from 
20.6 mag arcsec"^ to 22.6 mag arcsec"^. This compares to our 
20.5 mag arcsec"'^ and around 21.75 mag arcsec"'^ over the 
same range. 
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iDriver et all j2005h study the BBD in the MGC, using an im- 
proved SB definition involving elliptical isophotes and the deeper 
^app limit mentioned of 26 mag arcsec^'^ in B. Due to the overes- 
timation in the SDSS mentioned above, they find a syst ematic off- 
set of 0.4 mag to fainter values of the SB compared to lShen et alJ 
i2003h . They also see the mo re consta nt SB values at bright Mr but 
are otherwise consistent with lBlanto i et al. ( 2005c). 

An improvement to the study here would be to correct for 
the effects of internal extinction due to dust. This is particu- 
larly strong in highly inclined spiral discs. In this paper we at- 
tempted to reduce this effect by excluding galaxies with an axis 
ratio greater than an E7 elliptical using the axis ratio cut described 
in i|2.1l A rec ent study which does apply a dust correction is 
Ide Jong & Lacev, (.2000.). They u se a more local sample o f galax - 
ies from lMathewson et alJ h992h and lMathewson & Fori jl99d) . 
consisting of approximately 1,000 spirals of types Sb-Sdm. A de- 
tailed set of criteria that might be followed by future analysis are 
given bv lBovce & Phillippjl {\99^ . 

4.5 Colours 

The colours are computed using the Sersic magnitudes K-corrected 
to z = 0.1, the approximate mean of the sample. The r band shifted 
to 2 = 0.1 has Acff ~ 5600 A, similar to the V band. They are 
preferable to the observed frame colours as the galaxies over the 
whole redshift range are then more directly comparable. The LF 
fit implicitly assumes no evolution in the sample. One could take 
this into account by dividing the sample into redshift bins. All ten 
colours from ugriz were investigated and Figs. |6}|9| show the LF 
bivariate with u~r,g~r,r~i and r — z, a representative sample 
of the results. 

For the colours, the completeness (not shown) varies between 
60% and 70% for 13 < r < 16, dropping to 50% by r = 17.77 at 
the faint end and to consistent with zero at r ~ 11.5. 

In the LF split by colour, the expected trends of a fainter M* 
and steeper a for bluer colours are seen. There is a clear lack of 
faint galaxies with redder colours in all four plots. The panels split 
by colour show that the LFs are Schechter-like, with a steepening 
faint-end slope with bluer colour and possible upturns fainter than 
Mr — 51og/!. = —18, apart from in the bluest bins where the slope 
is already steep. 

In II — r we see a constant peak around u — r — 1.5 for 
-18 < Mr - 51og/i < -15. This then brightens to it - r ^ 1.75 
by Mr — 51og7i = —21 and disappears at brighter magnitudes. A 
red peak around u — r = 2.5 becomes evident at about the same 
magnitude and brightens to it — r ~ 2.75 by A/, — 51ogft = —22.5. 

The bimodality is clearer 'va.g~r, with the blue peak showing 
similar behaviour to that in it — r, moving from g ~ r ~ 0.35 to 
g - r 0.6 from Mr - 51og/t = -15 to A/,. - 51og/t = -21. 
The red peak is visible much fainter than in it — r, moving from 
g - r ^ 0.8 through 0.9 to g - r ~ 1 by A/,. - 51og/t = -24. 

In r — i, bimodality is no longer evident, with a single peak 
monotonically moving from r — i ~ 0.3 to r — i ~ 0.4 from 
Mr — 51og/i = —15 to Mr — 51og/i = —21. The peak is then ap- 
proximately constant at r — i ~ 0.4 at brighter magnitudes. There 
is therefore either no bimodality in this colour, or the peaks due to 
the populations causing the peaks in it — r and p — r are indis- 
tinguishable, either intrinsically or at the colour resolution probed. 
The constant colour brighter than Mr — 51og/i = —19.5 suggests 
that this could still correspond to the early type bulge-dominated 
population. The lack of an obvious blue peak is consistent with 
those peaks from it — r and g — r being due to star formation. 
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Figure 5. LF bivariate with absolute effective Petrosian 50% light radius 
surface brightness. The Petrosian radii are not seeing corrected. The distri- 
bution lacks bimodality and confirms the trend of fainter galaxies having 
lower surface brightness. The dataset contains an explicit surface bright- 
ness cut at ^RrQ = 24.5. The axis limits for the top panel are extended 
and the thick contours show the regions inside which volumes of 10^- 
lO'^ Mpc^ were probed. 



r — z shows similar behaviour to r — i, with the peak moving 
from r — 2 ~ 0.1 to r — 2 ~ 0.3, again constant brighter than 
Mr - 51og/i = -19.5. 

The plots for u — g, u — i and it — 2 are qualitatively similar 
to It — r, those for g — i and g — z to g — r and i — 2 to r — 
2, which suggests no obvious physical process apparent in these 
colours which has been missed by the four shown. 

The bivariate luminosity-colour distribution is also presented 
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bv lBlanton et aljj200lh and lOriver et alJ<200d) . lBlanton et all sub- 
divide the SDSS EDR by Petrosian g* — r* and see red and blue 
peaks at g* — r* ^ 0.3 and g* — r* ^ 0.75 respectively, with the 
same trends in the LF by colour. Our blue peak is the same but the 
red peak is at q — r ^ 0. 9 in both Petrosian and Sersic magnitudes. 

iDriver et alj i2006h show the u — r colour, divided at u — 
r = 2.1. They show the Petrosian value and that from the 
DRl PSF magnitudes, which are superior at separating the bi- 
modal population. For the core colours the Schechter parameters 
are M* = -19.15,-19.25, 0* = 0.0111,0.0136 and a = 
—0.15, —1.28 for the two populations respectively, thus showing 
a clear downward-turning but not Gaussian LF for the red galaxies. 

The prominence of the bluer population in u and g suggests 
that this peak in the LF is assoc iated with star formation. This is 
suppor ted by the GALEX LFs oflwvder et alj j2005h . lTrever et alJ 
J20()5h and lBudavari et alj<2005h . 

The u — r co lour is used in o t her in vestigations such as 
iBalosh et al] jgOO?) and 'Bal drv et all J2004h where the galaxies 
are clearly divided into two populations using counts in colour and 
morphology bins, which are fit by Gaussians. This is also done in 
B06 where u — r is compared to the neural network morphological 
type. 

Overall, the LFs subdivided by colour are consistent with a 
bimodal population of red bulges and star forming discs. 

4.6 Absolute Petrosian 90% Radius 

The bivariate LF is shown in Fig. 1101 The completeness is 50-60% 
for 13 < r < 16, zero at r 11 and 40% at r = 17.77. 

The expected pattern of fainter magnitude and higher normal- 
isation due to larger numbers of galaxies with decreasing radius is 
seen. There is a clear cutoff with magnitude for a given radius, giv- 
ing a declining faint-end slope at large galaxy sizes, which turns 
into a rising faint-end slope at small sizes. This suggests that the 
rising faint-ends in the Schechter LFs seen in other plots here are 
due to the smaller galaxies in a sample. As with surface brightness, 
the regions in which very small volumes are probed do not cut off 
the LF contours. 



4.7 eClass Spectral Type 

Fig. II II shows the bivariate luminosity-eClass distribution. The 
eClass is not measured in th e absence of a sp ectrum so there is no 
completeness map (although liall et alj \2Q04\ show that the eClass 
can be predicted to a RMS accuracy of ±0.06 using ANNs). The 
range —0.6 < eClass < 1.0 corresponds to early-late type galax- 
ies. 

Subdivided by eClass, the LF shows the expected bright- 
shallow, dim-steep faint-end slope corresponding to early-late type 
galaxies. For —0.2 < eClass < 0.0 the faint end LF is decreasing 
but not towards being consistent with a Gaussian LF. 

Subdivided by Mr, there is a clear peak at eClass ~ —0.15 
for galaxies brighter than Mr — 51og/i = — 18 which becomes less 
distinct at fainter magnitudes. There is a second much broader peak 
for Mr — 51og/i = —19.5 and fainter, moving from eClass ~ 0.15 
to eClass 0.25. 

The two peaks are thought, as with the colours above, to cor- 
respond to the passive bulges of bright ellipticals and star-forming 
discs. 

The LF has not been presented bivariate with eClass previ- 
ously in the literature. However, two recent studies use the similar 
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Figure 6. LF bivariate with u — r colour, extinction and K-corrected to 
z = 0.1. The expected trends of a fainter M* and steeper a for bluer 
colours are seen. 

PCA-based 2dFGRS rj class. This covers the range —5 ^ ^ 10, 
with low-high values corresponding to early-late types, in the same 
way as the eClass. Madawick et al. (2002) show LFs for 77 < —1.4, 
-1.4 ^ T] < 1.1, 1.1 < 77 < 3.5 and ri ^ 3.5. For the first bin, 
the LF shows a pronounced increase in normalisation fainter than 
Mr — 51og/i = —16, and appears Schechter-like in the other bins, 
with the usual dimming and steepening of the faint-end slope, al- 
though the Schechter function is always a poor fit. The 77 < —1.4 
clas s is better fit by a S chechter -1- power law function. 

iDriver et alJ llOOft) shows the LF in the MGC for the classes 
V = 1, r? = 2 and rj = 3,4, showing a similar pattern to 
iMadgwick et alj but without the faint-end upturns as the LFs are 
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Figure 7. As Fig.l6lbut for g — r colour. The results ai'e qualitatively similar 
to M — r with the red galaxy peak gaining strength. 
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shown to Mr — 51og/i < —16. Their MGC continuum type sepa- 
rates theirlJsjnore clearly. 

iKochanek et alj i200(i) discusses possible biases in spectral 
LFs due to the use of spectroscopic fibres for the latter. The spec- 
troscopic fibres will result in earlier spectral types for large galaxies 
due to bias toward their central regions. This bias is not quantified 
here, but their fig. 4 suggests that it is likely to be present even 
though most of our galaxies are poorly resolved ( i)2.1> . However, 
it is unlikely to change the broad behaviour of the LF split by Mr, 
i.e. a roughly constant low eClass peak and a higher peak which 
increases in value for fainter galaxies. Also, various studies (e.g. 
iGlazebrook et all2003h find that the effect is not large for star for- 
mation rate, which is a dominant constituent of the eClass and the 
?7 type via emission lines. 



4.8 Other measures 

Two further parameters were investigated: the galaxy stellar mass 
and the E jv surface density of galaxies. These are shown in Figs. 1121 
and ll3l again with no completeness maps as these measures require 
spectra. 

The stellar mass is useful as a physical parameter which could 
be directly compared with simulations. A tight correlation is seen 
between the LF and the stellar mass, with the LF in each bin show- 
ing a symmetrical shape of similar Gaussian form to that in absolute 
effective surface brightness. However, these results are expected 
because the luminosity is used in the calculation of the stellar mass 
in the first place, as part of the mass to light ratio calculated from 
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the models of star formation history (see 'Kau ffmann et alj|2003|l . 
For consistency with the other LFs, the r band is shown, although 
the correlation would be expected to be slightly tighter in i and z, 
due to these being less dominated by short-lived bright stars than 
the bluer bands. In the LP subdivided by stellar mass, the galaxies 
with a stellar mass of A/stdiar < 10^ Mq have distributions that 
are clearly truncated by the faint cut-off in absolute magnitude. This 
is evidence that there are further less massive galaxies not seen in 
this study. 

The S jv surface density of galaxies from the VAC is used in 
the morphology-density relation in B06. The bivariate LF appears 
to be of very similar form across all surface densities, with a power- 
law form, differing only in normalisation, which increases for both 
fainter magnitudes and lower densities. The faint luminosity cutoff 
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Figure 10. LF bivariate with absolute Petrosian 90% radius in kpc. The LF 
has a steeply declining faint-end slope at large galaxy sizes, turning into a 
rising faint-end slope at small sizes, suggesting that the steeper Schechter 
faint-end slopes seen in other plots are due to the small galaxies in the 
respective samples. As in Fig. |5] the axis limits in the top panel are ex- 
tended and the thick contours show the regions inside which volumes of 
10^-10^ h-^ Mpc^ were probed. 



due to the z > 0.053 limit on the density in the VAC is Mr — 
51og/i = —19.63, the top edge of the figure. The z < 0.093 limit 
does not have a strong effect on the bright luminosity galaxies. 

The variatio n in the LF w ith density and galaxy spectral type 
is investigated bv lCroton et alj ^05), who find that the LF is fit- 
ted by a Schechter function over all densities, that the (p* and M* 
vary smoothly over the density range but that the faint-end slope 
does not change much. Their normalisation takes into account the 
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Figure 11. LF bivariate with eClass PCA spectral type. There is a clear 
population of bright galaxies with low eClass which probably correspond 
to bright ellipticals. Some bimodality is evident. 
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Figure 12. LF bivariate with galaxy stellar mass. Although the z band abso- 
lute magnitude is likely to coiTelate with stellar mass slightly more tightly 
than r, the r band is shown for consistency with the other LFs. The correla- 
tion is expected as the galaxy luminosity is used in the estimation of stellar 
mass. 



volume in which galaxies are present in the density bin used. Our 
normalisation here is just the overall value for the bivariate LF, so 
the relative normalisations of the bins are essentially arbitrary. This 
is also seen in the LF divided by Mr, in which one expects fewer 
galaxies at low density with normalisation weighted by volume. 

It is also possible that the broad lack of change in the 
shape of the LF with density is related to the similar lack of 
change in slope in the correlation fun ction divided by l u minos - 
ity seen by e.g. IConnollv et alj i2002h and IZehavi et alj <2002h . 
but this would requ ire further investigation, especially now that 
IZehavi etal]j2004j) describes departures from the power law. 



5 FUNCTION FITTING TO THE BIVARIATE LF 

The Schechter function (equation |2j can, by taking L/L* 



(t)(M) dM = 0.4 In(lO) (f)* 



X exp [--10 



p.4(M'-J\/)i 



-10' 



0.4(M' -J\/)iQ + l 



dM. (24) 



If one then applies Bayes' theorem to the bivariate function then it 
can be written as 

7/i(M, X) dM dX = (t>{M) 'p(X\M) dM dX, (25) 
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Figure 13. LF bivariate with S5 surface density of galaxies. Unlike the 
other bivariate LFs here, the normalisation is arbitrary as it depends on the 
binning of the galaxies in density. 



which gives two-dimensional functions to fit to the SWML esti- 
mate. In general any well-behaved function can be fitted, expressed 
using either side of equation |25l 

One choice is to use a Schechter function for (j>{M) and a 
Gaussian with mean A{M — M*) + B and standard deviation ax, 
where A and B are c onstants for ip{X\M). This was first used by 
ICholoniewsk3jl985l) and is known as the Choloniewski function: 



ip{X\M) dM dX = 



X exp 



{X - A{M 



M* 



Bf 



2o\ 



dM dX. (26) 



This gives a 6-parameter fit. lCholoniewskil)l985l> actually used M 
rather than M — M* but the latter scales better with our data. 

Due to this greater dimensionality than the monovariate LF, 
a simple maximum likelihood fit is impractical, as the parameter 
space to explore is large. The method adopted here is the sim- 
plex search, which, while not of the greatest accuracy near the 
minimum, is robust to starting conditions and finds the minimum 
rapidly. 

Here the Matlab implementation, known as fminsearch, 
is used. It impleme nts the Nelder-Mea d simplex algorithm 
jNelder & Meadi l965h using the method of lLagarias et alj j 19981) 
as described at http://www.mathworks.com Options used were the 
Matlab default tolerance of 10~* and a maximum number of it- 
erations of 10000. The search can sometimes become stuck, for 
example, if the simplex decays to a lower number of dimensions 
due to the new point being too close to the centroid. However if the 
search converged it did so in less than this number of iterations. 

The measure of goodness of fit, the error function, is the stan- 
dard x^, given by 



X 



Ipobs — ^/i 

T err 



(27) 



The fit is done to log^Qi/) rather than iji so the error is estimated 
using 



Verr)] , (28) 



with i/Jerr from the information matrix determined as described 
above. 



6 FUNCTION FITTING RESULTS 

It was generally found that the functions are a poor fit to the data, as 
the deviation of the data from a smooth function is large compared 
to the size of the error bars, giving a large x^- Fig. 1 141 shows the 
fit of the Choloniewski function to the BBD and Fig. llSl the 
values. The fit is particularly poor for bright galaxies around Mr — 
51og/i ~ — 2L25 and /iRgo ~ 19.25 as the distribution narrows. 
The best-fitting parameters are cj>* = 0.0018, M* — —20.6, a = 
-1.19j£ji^^= OTL A = 0.30, B = 26.5. 

iDriver et alJ i2005h also find a poor fit, but for faint galaxies. 
They suggest that this is due to the assumption of a constant width 
for the Gaussian in absolute effective surface brightness whereas in 
fact the distribution broadens at fainter values. The assumption of a 
linear relation between the mean of the Gaussian and A/ may also 
be a factor. It is unlikely that the simplex is always getting stuck in a 
local minimum far from the global minimum as the fit is poor from 
a variety of starting points but the minimum reached is the same. 
One could add further parameters but six is already a fair number. 
The results here are consistent with the constant Gaussian width 
approximation also being the problem. We confirm the suggestion 
of Driver et al. ( 2005) that bulge to disc decomposition is needed. 

It is possible that no simple function will provide a good fit, 
as the superposition of underlying processes generating the bivari- 
ate LF may be too complex at the level of precision now available 
from the data. Or it may be the case that a single function will not 
fit, but the sum of two or more, each of which describes a different 
population, for example the red and blue or bulge and disc indicated 
above, may provide a better fit. The population might be subdivided 
usin g an optima l colou r- separation criterion, for example that used 
bv lBaldrv et alj i2004) . then fitting, for example, a Choloniewski 
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4>{x) = exp {~xy2)/V2n ^{ax) = 



(l){t)dt. (30) 



Figure 14. SWML estimate in greyscale and fitted Choioniewski function 
as contours for absolute effective Petrosian surface brightness. The right 
hand panel gives the LF logio ip value represented by the greyscale. The 
contours were fitted over the bins which contained five or more galaxies 
to minimize the effects of shot noise. The function provides a poor fit for 
bright galaxies. 




- 5logh [mag) 

Figure 15. values showing that the fit is generally poor, particularly for 
galaxies around Mr — 51og/i = —22 and fJ-R^g = 20. The right hand 
panel shows the values represented by the contours. 



function to each half. The difficulty is that the two distributions 
overlap, which would distort the two halves, and the sum of two 
Choloniewski functions involves a large number of parameters, es- 
pecially if the relation is non-linear between the mean of the Gaus- 
sian and AI, the Gaussians vary in width, or skew-Gaussians are 
required. 

Other bivariate LFs such as that with eClass were investigated 
with Schechter-Gaussian and Gaussian-Gaussian for earlier SDSS 
dat releases but the fit was also found to be poor. A skew-Gaussian 
was also tried using an available Matlab port of the software de- 
scribed at |http://azzaliru.stat.unipd.it/SN| The skew-Gaussian is 
given by 

fix) = 2(l>{x)^ax), (29) 
where (j>{x) is the ordinary Gaussian and ^{ax) adds the skewness 



A skewness of a = corresponds to the ordinary Gaussian, and a 
negative a gives the same shape as the positive but reflected about 
the vertical at a; = 0. The fits obtained were very similar to those 
for the ordinary Gaussian, with the skewness tending to low values 
of order 0.01. 



7 DISCUSSION 

A major question raised by these results and other bivariate LF 
studies is whether a theory or simulation can be constructed which 
will correctly predict the distributions over the different second LF 
parameters presented. Two major themes in the LFs are the bi- 
modality of the galaxy population and whether the underlying LFs 
per galaxy type are Schechter, Gaussian or some other form. 

Here we fitted a Choloniewski function to the bivariate 
luminosity-surface brightness distribution, and for earlier SDSS 
data releases experimented with it on other bivariate distributions 
and fitted Schechter, Gaussian, dual power-law Schechter and a 
sum of two Schechters to some of the individual LF slices. In all 
cases, the formal fits were found to be poor, with X^/^^ values of 
10 or more. It is thus difficult to tell when a fit is the best that one 
can do and quote the numbers or if the fit really is poor and should 
be improved. We therefore have not quoted function fit parameters 
for individual LF slices in this paper, although with more work one 
could attempt to do so. 

Concerning a bimodal population, the bivariate distributions 
are consistent with an early type, bright, concentrated, red popu- 
lation and a late type, faint, less concentrated, blue, star forming 
population. This idea of bimodality has been explored by others 
as it suggests two major underlying physical processes. These pro- 
cesses are thought to be connected to the formation of the bulge 
of the galaxy via mergers and the disc by accretion. Thus simula- 
tions of galaxy formation could be compared to the bivariate LFs 
presented. This is discussed further below. 

Several results support the idea that the underlying LFs of the 
giant galaxies, i.e. the usual E, SO and spirals are Gaussian, whereas 
those f or th e dwarfs are Schecht er. Sandage, Bin ggeli & Tammand 
il985l) and 'jerien & Tammann! Il997i) study the cores of nearby 
clusters;|de LaDoarent (2003) re views the LFs subdivided by mor- 
phology in various surveys and de Lapparent et alJi2003l) focusses 
on the ES O-Sculptor survey. A phy sical model supporting the idea 
IS given m ISchaeffer&sii9<1988l) in which the maximum galaxy 
mass is limited by baryonic cooling within a dynamical time-scale 
and a minimum mass is given by the minimum virial temperature 
required for cooling to occur. This gives the Gaussian LF for giants 
and gas stripping then g ives the Schechter function for dwarfs. 

iDriver et alJ i200d) find that the morphological classes E/SO, 
Sabc and Sd/Irr have distinct distributions relative to the bimodal 
population defined in the space of core it — r colour and Sersic in- 
dex n. The E/SO and Sd/lrr correspond to the red and blue peaks 
respectively, whereas the Sabc spread across the two peaks. This is 
argued to be strong evidence that the bulges and discs are the funda- 
mental two components to the bimodal population, not two separate 
galaxy types. Our ANN types contain few Sd/Irr due to the lack of 
these types in the training set, but our morphological LF supports 
this idea in the sense that the LFs for T types 1.5 through 3.5 (Sa 
through Sc corresponds to types 2-4) are essentially indistinguish- 
able, in contrast to the E/SO LFs which are distinctly brighter. 
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The bimodality in colour compared to that in morphology is 
also considered by B06 in the context of galaxy colour and ANN 
T-type morphology. There we find that whilst the distributions in 
colour are clearly bimodal and well-fit by a sum of two Gaussians, 
those in morphology are less so, with galaxies bridging the gap. 
This supports the non-fundamental nature of the Sabc types in this 
context. 

Also discussed by, most recently. Driv er et al.l \2Q0&) and 
lllbert et al. 1 .2006.) are blue spheroidal galaxies. These mimic faint 
ellipticals in morphology but are indistinguishable from blue galax- 
ies by spectra. These galaxies could cause type contamination in the 
LFs, for example causing a st eeper faint-e nd slope in the LF sub- 
divided by u — r, as shown bv Driver et al .. Again this agrees with 
the trends seen in our LFs and therefore requires spectroscopic and 
morphological information to resolve. 

Numerous studies ("e.g. IPopesso et af]|2005t fOonzalez et alj 
[2OO6) find that cluster LFs are better fit by a sum of a Gaus- 
sian at bright magnitudes and a Schechter at faint magnitudes, 
rather than a single Schechter function. The functions have a dip 
at Mr — 51og/i ~ —17, caused by the flattening of the Gaussian 
before the power-law faint end kicks in. Field galaxies are consis- 
tent with a Schechter function. A hint of a similar dip is seen here, 
as might be expected as the sample encompasses all galaxy envi- 
ronments. The faint end slope is thought to be caused by numerous 
faint red galaxies which are present in clusters but not in the field. 
Our LFs do not probe deep enough to say anything definitive. 

For the more Schechter-like results such as ours, the ques- 
tion is whether the LFs are as they appear, or whether the under- 
lying types are Gaussian, appearing Schechter because type mix- 
ing makes the functions tend towards the overall LF, which is 
Schechter, or a superposition of Gaussian at the bright end and 
Schechter or power law at the faint end. The idea has previously 
been discussed in terms of the Hubble morphological types and the 
effects can be seen here, particularly in Fig.|2| The LF for the early 
types could either be Schechter or Gaussian with a faint-end slope 
contaminated by dwarfs . In their sample which is essentially a sub- 
set of the one used here. lNakamura et alJi2003l) use the concentra- 
tion index to show that the sample is not contaminated with faint 
galaxies with soft cores, which will be dwarfs. Their concentration 
index LF is similar to their Hubble type LF. The same is seen here. 

The bimodality supports the idea that concentrated, early spec- 
tral type, red, high Sersic index galaxies correspond to the early 
morphological types. If this is the case then the Gaussian LF idea 
is also consistent with the LF plots bivariate with concentration in- 
dex, eClass, reference frame colours and Sersic index. These all 
show LFs for the early types which could be Gaussian with a rise at 
faint magnitudes, rising to a Schechter function for later types. The 
idea is also supported by the LF bivariate with absolute Petrosian 
radius, which shows steeply declining faint-end slopes which be- 
come less steep for galaxies smaller than ~ 2.5 kpc. 

Whatever processes cause the bimodality do not, however, 
show up in the LF subdivided by absolute effective surface bright- 
ness. However, it may well be the case that the sum of two bivari- 
ate distributions is still required and that, unlike the Choloniewski 
function, the relation between the mean of the distribution on the 
surface brightness axis must be a non-linear function of absolute 
magnitude to fit the lack of variation in the peak surface bright- 
ness at bright magnitudes followed by dimming at fainter levels. 
At fainter levels the two distributions could diverge and thus repro- 
duce the broadening seen. It may even be the case though that to 
get a formally good fit in the sense requires the distribution to 
be able to change width in the surface brightness axis as a function 



of absolute magnitude. If the Schechter function is also required 
to be a dual power law or Schechter plus Gaussian the number 
of parameters to fit rises to around 20. This may be impractical, 
at which point one would want to consider summing perhaps a 
greater number of simpler physically motivated functions, becom- 
ing more akin to semi-analytic models. These functions might be 
based on the Kormend y relation for bulges jKormendv 1977), Free- 
man's law for discs ( Freeman or the more recent relation of 
de Jong & L acey (2000) . and a relation for dwarf spheroids such 
as that in lMated Jl998l) . Another possibility is that central surface 
brightness may better distinguish between spirals and ellipticals 
than the absolute effective surface brightness used here. 

A useful way to connect observations to simulations is the 
conditional luminosity functi on, i.e. the luminosity distributi on of 
galaxies in a dark matter halo. lCoorav & Milosavlievi3i2005h con- 
struct the Schechter form of the galaxy LF using a model in which 
the central massive galaxies in haloes are distributed with a log- 
normal scatter and the satellite galaxies, where the total galaxy lu- 
minosity is greater than the central galaxy, are a power law. This 
predicts that the faint end of the LF is determined by the galax- 
ies in low mass haloes and should have a slope independent of 
colour selection or band. If this was combined with a prescription 
for predicting surface brightness then one might be able to predict 
the bivariate LFs presented here. This could be done by relating 
the absolute half-light r adius to the d imensionless spin parameter 
of the dark matter halo (Peebles 1969), as done by several studies 
(e.g. Fall & Efstathiou 1980; Dalcanton et al. 1997). If the same 
dark matter haloes are used for both models then one could pre- 
dict the bivariate brightness distribution, at least for discs. Relations 
of the haloes to other galaxy properties could simila rly predict th e 
other bivariate LFs, perhaps along the lines of e.g. lMo et alJi2004l) . 
However, our results show an LF faint end that does depend on the 
colour, so further comparison of observations and theory is clearly 
important. 

The conclusion from various studies is that the two physical 
processes are likely to be merger to form spheroids and accretion 
on to discs to give the star formation. In the sense that we see the bi- 
modal early-red versus late-blue population, our results support this 
conclusion. To test the ideas further requires bulge to disc decom- 
position and quantitative comparison with simulations incorporat- 
ing the ideas discussed above, as well as the future work detailed 
below. 

However, there is likely to be more to the story than 
simply mergers and accretion. For example, pseudobulges (e.g. 
.Korrnendv & Kennicu tt 2004) resemble bulges in disc galaxies but 
are formed due to internal (secular) evolutionary processes. Thus 
the relative importance of secular evolution on the observables dis- 
cussed here is a further avenue of investigation. Also, for example, 
none of the morphologies here take into account the presence of 
bars, which are important in a secular evolution s ense, if not in 
larger scale environment (e.g. Ivan den BerglJl2002h . There is then 
the question of more complex light profiles such as the core-Sersic 
model of Graham etal. (2003), types that don't fit the classifica- 
tions used here such as cD galaxies, higher redshifts, and so on. 

Important foUowup work to the work presented in this paper is 
to more rigourously separate the dwarf and giant galaxies in the LF, 
to apply reliable bulge to disc decomposition and to fit functional 
forms to the two underlying distributions. More refined structural 
parameters than those used to generate the Main Galaxy Sample 
such as fully Sersic-based measurements using elliptical isophotes 
should be used. Here dust-obscuration is addressed by simply re- 
quiring that the axis ratio is less than that for an E7 elliptical galaxy. 
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therefore excluding edge-on spiral discs that are highly reddened by 
internal extinction. Multiwav elength data could im prove this, for 
example matches to SWIRE I Lonsdale et al. 2003) data would be 
less affected. One could also match to GALEX ( Martin et al..2005i) 
data to extend to the UV for the star formation. 



8 CONCLUSIONS 

The conclusions of this paper are as follows. 

The luminosity function is measured bivariate with various 
galaxy properties from the SDSS and its associated value-added 
catalogues. The properties are eyeball morphological type, neural 
network morphological type, inverse concentration index C/inv, 
the Sersic index of the light profile, absolute effective surface 
brightness, the reference frame colours u — r, g — r, r — i and 
r — z, the absolute Petrosian 90% radius of the galaxy, eClass PCA 
spectral type, the stellar mass and the E5 surface density of galax- 
ies. Some of the parameters plotted here are new in the context of 
a bivariate LF, and for those that have previously been studied, the 
trends seen are confirmed at high signal-to-noise. This is the largest 
number of parameters studied in one framework for a bivariate LF. 

The SDSS is one of the best datasets currently available for 
this work, and the large sample combined with CCD photometry 
is ideal for precision bivariate LFs, enabling division of the galaxy 
population into numerous smaller samples, each of which is still 
highly statistically significant. In addition, the stepwise maximum 
likelihood estimator used is well-known and has been shown to be a 
good estimator. The magnitudes are K-corrected for each individual 
ga laxy using the weU-te sted K-correct code written for the SDSS 
bv lBlanton et"d] <2003th . 

The K-corrections do not take into account galaxy evolution or 
the effects of reddening due to dust, however we restrict the sam- 
ple to galaxies al z < 0.15 and those with an isophotal axis ratio 
less than 10/3 to minimize these effects. Near-infrared datasets of 
comparable size to the SDSS should improve the situation. 

The variation of the LF with all of the parameters is highly 
significant, both in shape and normalisation, except for S5 where 
the LF does not change significantly in shape and the normalisation 
is arbitrary. 

The bivariate LFs are consistent with an early type, bright, 
concentrated, red population and a late type, faint, less concen- 
trated, blue, star forming population. This idea of bimodality has 
been explored by others as it suggests two major underlying phys- 
ical processes, which in agreement with other studies we conclude 
are likely to be mergers leading to bulges and accretion on to discs. 
The bimodality is considered further in B06 in the context of galaxy 
colour and morphology as a function of environment. 

The various morphological measures are consistent, with the 
Sersic index showing the same broad trends as the C/inv, eyeball 
and neural network morphological type. 

The overall patterns of the LFs in the four colours u~r, g — r, 
r — i and r — z are similar, as are those for the other six from the 
SDSS ugriz bands. The main change is the increased star forma- 
tion in the bluer bands and the corresponding lack of bimodality in 
the redder bands. 

The LF bivariate with absolute galactic radius shows a steeply 
declining faint-end slope which becomes flat then rising for the 
smallest galaxies, supporting the idea that it is the small galaxies 
which cause the LF to be shaped like a Schechter function. 

The LF bivariate with galaxy stellar mass shows a tight cor- 



relation of higher mass for brighter magnitude, with Gaussian LFs, 
as expected as the luminosity is used in the calculation of the stellar 
mass. 

The LF bivariate with E5 surface density of galaxies shows 
little variation in shape over a large range in density from that in 
a cluster core to the field, being of power-law form throughout. 
Further work would be needed to quantify the variation that is seen. 

No definitive conclusion is reached as to whether the under- 
lying LFs for individual types are Gaussian or Schechter, due to 
the possibility of type mixing. This is particularly prominent in 
the LFs bivariate with eyeball morphological type assigned by the 
neural network, as those for the normal galaxies E, SO, and spirals 
might be expected to be Gaussian but if type mixing is present this 
causes them to tend towards the overall LF, which is well fit by a 
Schechter function. The idea is supported by the LFs bivariate with 
absolute galactic radius and to a lesser extent by Sersic index, con- 
centration index, spectral type and colours (the bimodality suggests 
that these correspond to similar galaxy populations). However, the 
plots are also consistent with Schechter functions with more or less 
steeply declining faint-end slopes. Further work would be required 
to investigate the effect of blue spheroidals, which mimic faint red 
galaxies in morphology but have bluer colours. 

It may be that the morphology of whole galaxies, particularly 
Sa-Sc spirals, is an intrinsically 'fuzzy' measure, so that the Gaus- 
sian LFs will not be seen. However, the contamination may be from 
dwarf galaxies so an important step is to separate these from the 
normal galaxies and determine their LF separately. The LF by ab- 
solute galactic radius is a step towards this but the next step should 
be using the azimuthally averaged light profiles available for all 
galaxies in the SDSS. Similarly, further work would be required 
to see to what extent the faint-end upturns seen in some of the LF 
slices are real. 

Some of the LFs are fitted with analytical functions, which 
are generally poor fits, probably due to this detail being beyond 
that reproducible by a simple function. An example is the fit of 
the Choloniewski function to the bivariate brightness distribution in 
which the fit may need to also allow for a variation in the width of 
the Gaussian with surface brightness, and a possibly non-linear re- 
lation between mean absolute magnitude and mean surface bright- 
ness, which results in the addition of two parameters to the six al- 
ready present. The poor fit is also seen in fDriveretalJj2005l) . where 
it is interpreted in terms of the diverging trends of bulges and discs 
with decreasing surface brightness. The fitting method used is eas- 
ily extended to any of the other bivariate LFs presented here. 
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APPENDIX A: COMPLETENESS MAPS 

For each result in which the second LF parameter does not require 
spectroscopy, i.e. Tjpg, Tann, C/inv, and n, the completeness 
maps for the bins used in the LF SWML estimate are shown. The 
parameters are therefore selected with the explicit cuts implied by 
the extent of their axes on these maps. For /^Rgo we show the com- 
pleteness as a function of apparent SB. Maps for the other param- 
eters are shown as a function of apparent magnitude. With the ex- 
ception of Tjpg, which is selected from a smaller area with much 
higher spectroscopic coverage, the overall completeness will gen- 
erally not exceed the ratio of the spectroscopic areal coverage to 
that of the imaging over the whole survey. This is about 70%, from 
4,783 or 4,681 deg^ersus 6,851 or 6,670 deg^ for the VAGC and 
the CAS respectively. In general, the completeness will be lower, 
due to further objects being masked out for the various reasons de- 
scribed above. The maps give an indication of the possible extent 
of biases in the LFs due to incompleteness in spectroscopy. 

Whilst these maps show the completeness in spectroscopy 
compared to selected targets, they do not show incompletenesses 
in the imaging compared to the true galaxy distribution. These will 
arise if either the object is not detected, or it is masked out from the 
photometry that is compared to the spectroscopy in the maps. 

The photometric incompleten ess due to su r face bri ghtness 
selection effects is investigated by iBlanton et alj i2005dl for the 
ranges in apparent surface brightness and magnitude 18 < /iapp < 
24.5 and 14 < r < 17.5. They process simulated galaxies with 
the same software used to process the real DR2 imaging data. The 
effects of seeing, noise and sky subtraction are included. At low 
surface brightness, they find that the overall completeness drops 
below 90% around /iapp > 22, reaching 50% at /^app = 23.4 
and at /^app = 24.3. This drop is fairly uniform over the whole 
range in r. At the bright end, there is slight incompleteness (90%) 
for dim high surface brightness objects (/iapp < 19 and r > 17), 
ascribed to the star/galaxy separator misclassifying the objects as 
stars in marginal seeing. There is more substantial incompleteness 
for bright high surface brightness objects, the overall completeness 
dropping to 80% for /iapp > 18.5 and 64% for /iapp ~ 18.2. 
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This is ascribed to the explicit surface brightness cut in the Main 
Galaxy Sample to eliminate binary stars f i|2.U . The overall frac- 
tion of galaxies affected by incompleteness at the faint end is small 
because the numbers observed drop well before the completeness 
drops. 

At the bright end, objects will be missed either through being 
shredded into smaller objects by the deblender or, as mentioned, 
by being misclassified as an object other than a galaxy. Besides the 
marginal seeing example mentioned, in recent years evidence has 
also emerged for a class o f galaxies, ultra-compact dwarfs (UCDs, 
e.g. IPhilUpps et aT]|200lh . which are affected by this. Again any 
effect on the overall LF for the apparent magnitude range investi- 
ga ted in this paper is un likely to be significant because, as shown 
bv lBlanton et aljfcoOSd) . few objects are observed in this regime. 

The VAGC quality flag rejects objects which are bad deblends 
in the ranges z < 0.01 and Mr — 51og/i > — 15. Most bad deblends 
are likely to be fragments of large galaxies and would therefore be 
within this range, being both nearby due to the apparent size of the 
galaxy in which they are contained and intrinsically faint due to 
their small mass and therefore light compared to a galaxy. Objects 
without a spectrum will not be in the range but will also not be 
sho wn in the LF. 

IStrauss et alj i2002h show that in the range — 23 < Happ < 
—24.5, 35% of the objects are bad deblends after applying their 
local versus global sky value cut. Therefore objects that are in this 
35% but not in the range rejected by the quality flag may remain 
in the LFs presented here. Accurately quantifying the number of 
objects for which this is the case would require detailed inspection 
of images in large numbers of LF bins, which is beyond the scope 
of this paper. Here we inspect the images for some of the outlying 
bins and plot the LF slices (but not the greyscale) only for bins 
containing more than 20 galaxies, for which any contamination is 
small. Bad deblends are therefore unlikely to cause any significant 
biases in the overall LFs. An exception is again made for Tjpg, 
where the cut is 5 galaxies per bin, because these were manually 
inspected for outliers by the JPG and so the shot noise becomes the 
limiting factor. 

Objects missed altogether by the survey will mean that we are 
not m easuring the true LF. However, studies (e.g. iHavward et all 
l2005h have shown that low-surface brightness galaxies do not 
contribute significantly to the luminosity density of the universe, 
particularly in the range of brightness we study here. Similarly, 
UCDs have absolute magnitudes of —13 < Mb ^ —11 
jPrinkwater et all2000l) and so would not enter into our LFs. 
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Completeness 




Figure Al. Completeness values for the bivariate luminosity-TjpG distribution. The top left panel shows the function in a similar manner to the upper panels 
of Figs. [Till 31 The top right and bottom left panels show the completeness as a function of apparent magnitude and Tjpq respectively. The bottom right panel 
shows the V/Vmax statistic. The selection boundaries are the edges of the plot on both axes. 
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Figure A2. As Fig.|M]but for Tann- 
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Figure A3. As Fig.lATIbut for Ch 
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Figure A4. As Fig. lAll but for Sersic index n. 
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Figure AS. As Fig. lAll but for /^r^q . 
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